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Abstract 

A damping strategy for blisks (integrally bladed disks) of turbomachinery involving a 
friction ring is investigated. These rings, located in grooves underside the wheel of the 
blisks, are held in contact by centrifugal loads and the energy is dissipated when relative 
motions between the ring and the disk occur. A representative lumped parameter model of 
the system is introduced and the steady-state nonlinear response is derived using a multi- 
harmonic balance method combined with an AFT procedure where the friction force is 
calculated in the time domain. Numerical simulations are presented for several damper 
characteristics and several excitation configurations. From these results, the performance of 
this damping strategy is discussed and some design guidelines are given. 

Keywords: Bladed-disk, friction damping, nonlinear dynamics, harmonic balance method 
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1 Introduction 

Turbomachinery bladed disks operate in severe environments in terms of aerodynamic loads and 
may experience high vibratory stresses due to resonance or flutter. Their ability to withstand 
these stress levels, which may cause high cycle fatigue, mainly depends on the external inputs 
of damping. In bladed disk assemblies, friction in interfaces is the most widely used source of 
external damping; this includes shroud contact in blades or platform dampers. These devices 
were widely studied in the past and numerous methods of design and analysis are found in 
the literature. The interest reader c an find examples based o n lumped param e ter or equ ivalent 
models in works of iGriffinI (Il990l b ISinha and GriffinI (|l983l ). IPerri and Hed3 (Il998l^ or ICsabal 



1998 ) and, more rece i itly, e xte nded to more comp l ex str uctures in works by iGuillen and Pierre 



19991 ). iNacivet et all (|2003l ) or IPetrov and EwinsI (|2003l ) among others. The contact interface 



could be on e dimensional, two diniensio nal or three dimensional (e.g. Sanliturk and EwinsI . 1996 : 
Ghen et al.l . bood : iNacivet et all . l2003l ) 



pesiae, numerous e x perimental works were presenteu botn m contact Rmematic description 
fseelMeng et al. ( 1986): |Ya,ng and Mena| (1998) for example) and dynani ic response prediction 
3.g.TBcrthi llier et al.l . Il998l : IVang and Menal . Il99i : ISanliturk et~all . boOlh . 



However, in blisks (integrally bladed disks where disk and blades are a single piece), inherent 
damping is very low and mechanical joints no longer exist leading to the disappearance of nearly 
every sources of energy dissipation. Consequently, new methods and devices need to be found 
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for the damping of blisks. Among others, friction damping using a ring is a solution which is 
found, experimentahy, to be efficient on single piece rotating structures (such as labyrinth seals 
for example). These rings are located in dedicated grooves underside wheels and contact is 
permanently maintained by centrifugal load due to engine rotation; friction and slipping occur 
during the differential motion of the two bodies in contact. 

As o pposed to under-platform dam pers in bladed disk assemblies, these devices were seldom 
studied; Ziegert and Niemotka ( 19931 ) have presented methods of analysis and optimal design 
of split ring dampers for seals were the motion is studied using a quasi-static beam-like de- 
scription (strength of materials). This quasi-static description lead to a good understanding of 
the phenomenon and give some qualitative results. However, a nonlinear dynamics analysis is 
required to raise satisfying design rules. 

In this paper, methods for predicting nonlinear steady-state response of blisks with friction 
ring dampers are presented along with some qualitative results based on a representative lumped 
parameter modelling. The dynamical analysis is performed using the multi-Harmonic Balance 
Method {HBM). This method is found to be efficient and was widely used, in the past, to study 
friction dam pers in bladed di s ks ass e mblies or shrou d cont acts. Some examples have been 
presented by Wang and Chen ( 1993 ). Sanliturk et al. ( 1997) and Yang and Menol (1997) for 
-harm onic vibrations and by Pierre et al. (jl985l ). Petrov and Ewins ( 2003 ) . iNacivet et al. 



mono- 



(|2003l ). and lGuillen and Pierrd (|l999l ) for multi-harmonic vibrations. 

Results from numerical simulations (using this nonlinear analysis) are then presented and 
described. Several parametric studies are presented in terms damper characteristics, contact 
characteristics, excitation configuration and mistuning are presented and some design guidelines 
are raised. 



2 Lumped-parameter model of a blisk with friction ring damper 



Lumped parameter models have been widely used to study friction. Beside from being sim- 
pler to study and compute (particularly for parametric studies), they can provide some inter- 
esting and valuable qualitative results if correctly designed. In turb omachinery applications , 
se veral models have been proposed. Among others, we can men tion Sinha and Griffin ([l98 j ) 



„_ Wang and ChenI hmi ) for friction damping applications and iGriffin and Sinhal (|198 
Lin and Mignoletl (|l996l )" for mistuning applications. 



and 



In this paper, we present a new lumped parameter model. It consists in an elementary sector 
of a blisk and its associated friction ring element and it is represented in figured! 

As this system is rotationally periodic (or cyclic), its modes are such that the deflection 
of a given elementary sector (described in figure [1]) is the same as, but with a constant phase 
difference from, the preceding (or following) sector. If N is the number of elementary sectors of 
the structure, the phase angle is lim/N where n is called the nodal diameter number. There are 
N possible values of n (from to — 1) but only N/2 + 1 values (from to N/2), if N is even 
or (N — l)/2 -|- 1 values (from to (A — 1) /2), if N is odd, represent different mode shapes; the 
others being orthogonal to them with the same eigenfrequency. Consequently, except for n = 
or n = N/2 if A^ is even, the natural freq uencies a. r e rep e ated. Further in formation on this 
subject can be found in the literature (e.g. iThomasl . Il979l : IWildheiml . Il98ll ). As an example, 
figure [2] displays the natural frequencies versus nodal diameter number of the (blisk) lumped 
parameter model with A = 24 sectors. 

The blisk lumped-parameter model was designed so that it stands as much as possible for 
a real blisk, particularly in terms of coupling between the blades and the disk. There are three 
families of modes and regions of various modal density. This modal density is directly linked to 
the coupling between the blades and the disk in the particular mode. A general rule, which for 
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Figure 1: Lumped-parameter model 
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Figure 2: Frequency / diameter diagram 
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that matter applies to the present example, is that the blade/disk coupling can be important is 
low frequency/low nodal diameter and high frequency/high nodal diameter regions whereas in 
low frequency/high nodal diameter or high frequency /low nodal diameter regions this coupling 
is generally weak. For example, in the two first families the coupling between the blades and the 
disk is weak for high nodal diameters whereas in small nodal diameters the coupling is stronger. 
Also note an important coupling appearing in the 2"'^ and 3'''^ modes with nodal diameter 
number 2. Recalling that, in operating conditions, the external forcing generally acts on the 
blades this coupling parameter is an important factor regarding the efficiency of the friction ring 
damping since it determines the disk participation in the global motion. The energy dissipation 
due to friction occurs when the relative displacement between the disk and the ring is enough; 
this is achieved when the blisk is excited in areas of strong coupling where both the blades and 
the disk move. The numerical values of the lumped-parameter are shown in table [TJ 



Table 1: Numerical values of the model 

mp ko kc rni kl ^ ^ 

1,2 kg 5.10'^ N/m 6.10^ N/m 350 g 1.10^ N/m 250 g 2.10^ N/m l%o 



Finally, concerning the ring, its mass and stiffness are chosen so that its lumped-parameter 
model represents a circular beam in extension motion and its stiffness k^ represents the longi- 
tudinal rigidity of the ring. 

An examination of the local behaviour of the friction area can help understanding the fric- 



tion dissipation mechanism. With reference to figure 3(a) , both the disk and the ring experience 
bending motion (bending stress at,) and are held in contact (slipping or not) due to the cen- 
trifugal load P. The friction forces in the interface (T < fiP) between the two elements will 
generate tension or compression stresses {(Tt/c) niainly in the ring which ensure the continuity 



of stresses at the interface as shown in figure 3(a) As a consequence the friction load is directly 



linked to the tension stresses and, assuming a linear elastic deformation, to the axial rigidity of 
the ring. 
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(a) Local behaviour in the frictional region 



(b) Friction ring in its groove 



Figure 3: Physical model and local behaviour 

Figure [3 (b) | shows a physical model of the friction ring in its groove. In order to be assembled 
in these groove, underside the wheel of the disk, the ring needs to be split; the influence of this 
split on the dynamics of the system will be discussed later. 

The system of equations of motion can be represented as: 

Mi(t) + Ci(t) + Kx(t) = fe^t(t) + ini{t) (1) 

where x = [x^, x^, x^^ , Xfjj]"^ is the physical displacements vector which gathers displacements 
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of the disk, the blades and the ring degrees-of-freedom. The matrices M, C and K are the 
structural matrices of mass, damping and stiffness respectively. The vector fg^i is an external 
distribution of rotating forces of Engine Order type, so that, considering a bladed disk with N 
sector, the j*'' sector's forcing term is: 



(2) 



where p is engine order. 

The vector of nonlinear forces f„i due to friction is expressed, according to the Coulomb law, 

as: 

ini{t) = -/xPsign(u(i)), (3) 

where u(t) = Xd{t) — Xr{t) is the relative displacements vector between the ring and the disk, /x 
is the constant friction coefficient, and P the normal contact load, which can be expressed for 
rotating structures as: 

P = rrirRn'^ (4) 

where rrir is the ring's mass, R is the distance from the axe of rotation and the ring and is 
the engine rotation speed. Due to the high contact load resulting from the engine rotation, is 
assumed that the ring and blisk remain in contact at all times. 



3 Nonlinear Analysis 

3.1 Harmonic Balance Method 

The aim of the nonlinear analysis performed here is to calculate the steady-state response of the 
system under a periodic external forcin g. Among the variety of method available to achieve thi s 
goal, the Harmonic Balance Method ( Navfeh and Mook . 1979 : Szempliiiska-Stupnicka . 199d ) 
presents several advantages in terms of efficiency and accuracy with respect to its computational 
costs. 

The idea is to express the time-dependent variables of the nonlinear equation ([T]) in terms 
of Fourier series (truncated in practice); the response x(t) is then: 



x(t) = X° + ^{X"'^cos(nwi) + X"'"sin(na;t)}. 



(5) 



n=l 



The new unknowns of the problem are the Fourier components X"''^ and X"'^ which are balanced 
(using a Galerkin procedure) and the equation of motion can be rewritten using this formalism 
for each terms of equation ([1]) in the frequency domain as: 



with: 
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The matrix A is block-diagonal, A = diag (K, Ai, . . . , Ajy^), and: 



-kujC 



kujC 

iujfm + K 



(8) 



The problem ([6]) is nonlinear and is usually solved iteratively using a Newton-like method. 

One of the issue of the HBM is that the Fourier components of the nonlinear forces F„/ can 
not, in most cases, be derived straightforward. In particular, in problems involving friction, 
the nonlinear forces are only known (using equation (l3|)) in the time domain and the expansion 
in the frequency domain is not obvious. When few harmonics (no more than 3 in practice) 
are retained an arialytical derivation of these frequency domain components is possible (e.g. 
Nayfeh and Mook , 1979 ; Wang and Chen , 19931 ) . To overcome this proble m when more harmon- 
ics ar e considered, the use of an Alternating Frequency Time method (e.g. ICameron and Griffinl . 
19891 ) can be a solution. This method will be detailed in the next section. 



3.2 AFT Method 

The aim is, in the Newton procedure, to derive the Fourier components of the nonlinear forces 
as a function of the displacements Fourier components. In computational applications, the idea 
of the AFT method is to use Discrete Fourier Transformation (DFT) to derive the Fourier 
components of the nonlinear forces for given displacements in the frequency domain. 

Let's start from X"'*^ and X"''^, the harmonic components of the response, predicted by a 
given iteration of the Newton-like method. Using, an IDFT {Inverse Discrete Fourier Trans- 
form) procedure, one can express the associated displacements x(t) and velocities x(t) in the 
time-domain. Then using a nonlinear operator the nonlinear force in time-domain is derived 
(here this step is achieved using equation ^ as described in the next section). Finally, using 
a DFT {Discrete Fourier Transform) algorithm one expresses the harmonic components of the 
nonlinear force in the frequency-domain. The following diagram illustrates this procedure: 



X"'' 



Xn,. x(t), ±{t) 



^ nl ' ^ nl 



DFT 



f„z(x,x,t) 



3.3 Calculation of friction forces in the time domain 

At a given iteration of the Newton scheme and during the AFT procedure, the friction forces 
need to calculated, in the time domain, in accordance with the velocity and displacement given 
(by the Newton method). This is generally done using the Coulomb law or, more usually, using 
a regularized Coulomb law and several strategies can be used. Some a uthors have proposed 
a regularized expression of the sign function (e.g. Berthillier et al. , 19981 ) or arctangent or hy- 
perbolic tangent f unction (e.g 



distributions fe.g. iPierre et al 



'etrov and Ewins . 20021 ). some others have used the theory of 



1985ll. and some have used an iterative procedure i n the time 



domain (e.g. Poudou and Pierre . 2003 : Guillen and Pierr3 . 1999 : Petrov and Ewins . I2OO3I ). For 
the latter strategy, it is convenient to add a tangent stiffness to each friction element; doing so, 
the Coulomb law can be used straightforward to compute at each time step the equilibrium of 
the contact point. In previous studies on under-platform dampers, this stiffness was included 
in the damper model. However, in general contact problem, this tangent stiffness appears as 
an equivalent penalty stiffness and can then account for the elasto-plastic shear deformations 
of the contact asperities of the bodies in contact. 

Figure m shows a friction element defined in relative displacements, with a tangent stiffness: 
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u is the relative between two points in contact, u = — Xj. in our example, 
z is the relative displacement of the contact point. 




Figure 4: Friction element (relative displacements) 
The nonlinear forces vector given by equation ([3]) becomes: 

f m = J-^t ("W - ^(*)) if ^tHt) - z{t)\ < fiP, 
"'^ ^ \-/iPsign(z(t)) if kt\u{t)-z{t)\>fiP. ^' 

In the time domain of the AFT procedure, velocities and displacements are given in one 
period of the motion. These are entries for the computation of the associated nonlinear forces 
in one period of the motion. This computation has to be performed time-iteratively (like a 
numerical integration) because of the history dependency of the friction law. For example, 
given the relative displacements and velocities at time ti, u(tj) and u{ti) and the displacement 
of the contact point z(tj), the friction force can be predicted (using equation ([9])) for time ti 
assuming a sticking state (z(ij) = z(fj_i)) for the contact point: 

f^l' = -kt{u{U)-z{U.i)). (10) 

The Coulomb law is then used to correct the friction force and the displacement of the contact 
point: 



nl I nl 



z{u) = uiu) - (lib) 



Using this iterative scheme, the periodic friction force can be derived. 



4 Numerical results 

In this section, several aspects of friction ring damping are investigated through numerical 
simulations and discussed. All nonlinear numerical simulations were performed using the multi- 
harmonic balance procedure and a third order Lagrange predictor was used for the continua tion 
of the frequency responses. A d*^ order Lagrange predictor, see Stoer and Bulirsch ( 1980l ) for 



example, requires the last d + 1 points to predict the next one X^"*"^ and an interpolation is 
build on the arc-length s (which is here the excitation frequency) as: 

. , , . A. .(p+i) _ .(p-d+j) 

= L, . X--- with L. = p (12) 

1=0 

The step As is adjusted a posteriori as a function of the (Newton-like) corrector's performance 
(number of iterations) at the previous step. 
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All frequency response plots represent the blade response (degree-of-freedom in figure 
[T]), with, in dashed line, the linear response (without ring) and in solid lines, the nonlinear 
responses (with friction ring). Parametric studies on the type of excitation (or blade/disk 
coupling as mentioned earlier), on contact parameters, damper characteristics are presented 
and the efficiency of this damping strategy is discussed. Then the influence of ring's split 
discussed. And finally, the impact mistuning and friction ring damping is investigated. 

4.1 Optimal damper mass — Blade/disk coupling 

In this first part, the influences of the ring damper mass and of the excitation configuration are 
investigated. The rotation speed is kept constant is these simulations; however, note that as 
the ring's mass changes, the contact normal load also change with respect to equation 

Strongly blade/disk coupling With reference to the frequency/nodal diameter map of 
figure [21 we first focus on a small engine order excitation (p = 2 for example in equation ([2])) 
with frequency range around the first modal family (similar behaviour would appears in other 
modal family however). 

Several frequency responses are plotted in figure [5] for different damper mass. From these 
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Figure 5: Frequency responses with strong blade/disk coupling; ( ) : linear response, ( ) : 

nonlinear responses at different damper mass (g) 

response, we can first notice that the damping is efficient in this configuration when the ring is 
sufficiently thick and that an optimal value appears {rriring ~ 20g). At this optimal value, the 
reduction of vibration amplitude is nearly 4. Also note the progressive frequency shift of the 
resonance peaks as the ring's mass increase, this is due to two facts combined. First, considering 
the associated linear problem, as the ring's thickness increases the global system becomes stiffer 
which account for the frequency shift to the higher frequency. And second, the nonlinear effect 



8 



of friction has usually a softening effect on system's dynamic and, with respect to figure O 
as the ring's mass decreases the contact normal load also decreases (assuming a fixed rotation 
speed) and more and more slipping occurs which account for the frequency shift to the lower 
frequency as the damper's mass decreases. 

Weakly blade/disk coupling A weaker blade/disk coupling can be found in higher nodal 
diameter in figure [2j As an example, let's consider a 6 engine order excitation on the second 
modal family and the associated frequency responses plotted in figure [6l In this example, the 
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Figure 6: Frequency responses with weak blade/disk coupling; ( ) : linear response, ( ) : 

nonlinear responses at different damper mass (g) 

friction ring damper is clearly less efficient; even if an optimal mass appears (around 5g), the 
reduction of vibration's level is not as important as in strong coupling configuration. In effect, 
in the configuration, there is not enough motion in the disk as the blades are excited and as a 
consequence no relative motion between the disk and the ring can occur. 

4.2 Friction ring damping efficiency at different operating speeds 

In this section, we focus on the infiuence of the normal contact load (which is, in our applica- 
tion, related to the engine rotation speed) on the friction ring damping efficiency. To do so, 
the damper mass is fixed as well as the excitation parameters (Engine Order 2) and several 
simulations were performed for different values of the engine rotation speed. 

In figure \7\ the results are depicted and the first remarks is that the existence of an optimal 
rotation speed/normal load is clear around 3000 rpm. Moreover, one can see that that in a 
quite large vicinity of this optimal value (from 2000 to 6000 rpm), the friction ring damping is 
still efficient. 

This is a key feature since a turbomachinery often operates at different rotation speeds 
(take-off, cruise, landing, ...). As a consequence, the rings can, for example, be optimally 
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Figure 7: Frequency responses at various rotation speed; (---): linear response, ( ) : 

nonlinear responses at different rotation speeds (rpm) 



designed for the cruise speed and be also efficient at take-off or landing speeds. 
4.3 Symmetry breaking induced by ring's split 

As mentioned earlier, in order to be assembled the ring need to be split. However, considering 
the model described in figure [H a split in the ring would necessary spread over one elementary 
sector and generate an exaggerated gap compared to the technological reality. As a consequence, 
we built a refined model, presented in figure [HI in which the disk mass, mo, is refined into three 
masses, TOq. This allows the split to spread only over a portion of one sector. 

The influence of the ring's split was then investigated by comparing the dynamic response of 
the system with or without a split. First, the ring is assumed to be fully stuck (no sliding occurs) 
and the system is then linear. The frequency response plots resulting from these simulations 
are displayed in figure O 

We focused on the first resonance peak of the 2 engine order excitation and comparing the 
blade's resonant responses with (figure |9(b)[ ) or without (figure 9(a) ) ring's split, it appears 



that each blades have a different response and that the maximum level is more important when 
the ring is split. This symmetry breaking, which corresponds to the separation of the initial 
twin modes (two-fold degen eracy), is similar to the consequence of mistuning and can lead to 
localization phenomena rsee lWei and Pierrj B H) for example). 



When the nonlinear effect are introduced (figure [TOl) , the same phenomena of symmetry 
breaking and increasing of resonant responses appear. However, the response amplification is 
limited by the nonlinear effect and then smaller than in the linear case. 

These results show the influence of the split in the ring damper on the global dynamics of 
a bladed disks. However, these observations have to be put in perspective with respect to the 
influence of the inherent mistuning (discussed next) of a real bladed disks which may generally 



10 




Figure 8: Fine lumped-parameter model 




Figure 9: Linear frequency responses with and without ring's split; ( ) : without a ring, 

( ) : with ring (split or not). 
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Figure 10: Nonlinear frequency responses with and without ring's spht; ( ) : without a ring, 

( ) : nonUnear responses. 

be more important. 

4.4 Influence of mistuning 

To conclude on the numerical results, the influence of mistuning was investigated. Mistun- 
ing refers to small variations in the properties of a system originally designed to be cyclically 
symmetric. These variations are mainly due to manufacturing tolerances or material inhomo- 
geneity and are inevitable. Mistuning can lea d to significant c hanges in dynamical behaviour 
with respect to the tuned case as presented by Wei and Pierre (jl988al lbl): vibratory amplitudes 



in forced response are increased and impact on high cycle fatigue can be considerable. These 
effects are even more important in Integrally Bladed Disks since their "natural" damping is 
far smaller than for bladed disks assemblies. The i nflu ence of mistuning on fric tion damping 
has been studied earlier by Griffin and Sinha ( 19851 ) or Lin and Mignolet ( 19961 ) and the per- 



formance of the damping technology discussed here as applied to blisks may be influenced by 
this phenomenon. 

Considering the lumped parameter model of figure [H mistuning was introduced as random 
perturbations in the blades' stiffness kl and k'^ as it is often the case since the disk is assumed 
to be symmetric. These perturbations lead to changes in the natural frequencies of the blades 
(cantilevered) . 

A frequency response is shown in figure \TT\ the mistuning standard deviation is about 1.2% 
of the cantilevered blade's nominal natural frequency. The resonant responses of the blades are 
displayed for the tuned and the mistuned case with or without ring. First, we can notice the 
typical feature of mistuned responses which is the break of symmetry which leads to a separation 
of the resonance peak (with or without rim damping). Then, the vibratory amplitudes are 
increased in the mistuned case with respect to the tuned case; however, the damping is still 
efficient in presence of mistuning. 
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Figure 11: Frequency responses with and without mistuning; ( ) : hnear responses, ( ) : 

nonhnear responses. 

4.5 Conclusions on the friction ring damping performance 

The results of these numerical simulations clearly highlight the possibilities and limits of the 
friction ring damping strategy. In essence, the dissipation mechanism of friction is efficient only 
if some relative motion between the two bodies (the disk and the ring) in contact occurs. As 
opposed to under-platform dampers, the rings are not directly in contact with the blades to 
which the external forcing applies. As a consequence, the coupling between the blades and the 
disk which conditions the energy transfer from the blades to the disk is an essential criterion. It 
has been demonstrated that when the coupling is weak (that is when little energy get transferred 
to the disk), the ring damping is less efficient. On the other hand, when this coupling is 
sufficiently strong the rings are efficient. Finally, it was shown that the contact parameters 
(normal loads or friction coefficient) are also significant in the friction damping performance. 
The qualitative results given by these simulations are quite similar to those observed with under- 
platform damper; the alternation between slipping and sticking states results in frequency shift 
which offsets the resonant frequency from its initial value (i.e. non damped), and an optimum 
can be reached for a given excitation configuration in the ring's mass. 

5 Conclusions 

A strategy of damping for Integrally Bladed Disks (blisks) using a friction split ring was inves- 
tigated. A procedure of analysis using a multi-Harmonic Balance Method and an ^FT method 
was described and applied for calculating the steady-state nonlinear response of a system under 
periodic excitation. 

Simulations on the lumped-parameter model presented lead to qualitative results in terms 
of energy dissipation with respect to several configurations. Some guidelines on the use of this 
technology, its limits and performances, were highlighted. We show that the coupling ratio 
between the disk and the blades was a critical parameter of the damping power of such rings. 
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